Spin-flip scattering and critical currents in ballistic half-metallic d-wave Josephson junctions 
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We analyze the dc Josephson effect in a ballistic superconductor/half-metal/superconductor junction by means 
of the Bogoliubov-de Gennes equations. We study the role of spin-active interfaces and compare how different 
superconductor symmetries, including d-wave pairing, affect the Josephson current. We analyze the critical 
current as a function of junction width, temperature, and spin-flip strength and direction. In particular, we 
demonstrate that the temperature-dependence of the supercurrent in the d^-symmetry case differs qualitatively 
from the s- and d x 2_ y 2 -symmetries. Moreover, we have derived a general analytical expression for the Andreev 
bound-state energies which shows how one can either induce 0— 7r-transitions or continuously change the ground 
state phase of the junction by controlling the magnetic misalignment at the interfaces. 
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I. INTRODUCTION 

The antagonistic nature of superconductivity and ferromag- 
netism makes their coexistence an unlikely one in bulk mate- 
rials. In conventional superconductors, the current is carried 
by Cooper pairs consisting of two electrons in a spin singlet 
state, i.e. with anti-parallel spins. 1 Ferromagnets on the other 
hand, favor parallel spin-alignment. Hence, one might expect 
that the proximity-effect arising when bringing a ferromagnet 
into contact with a superconductor would decay rapidly due 
to this antagonistic nature. 

However, it has been predicted that this rapid decay may 
not occur if there is some magnetic inhomogeneity present 
at the interface between a superconductor and a ferromagnet. 
Bergeret, Volkov and Efeto\ 2 proposed a theory that coherent 
triplet Cooper pairs can be induced in ferromagnets, which in 
turn gives rise to a long-range triplet current. Several experi- 
ments have detected signs of this long range triplet current, 3 7 
and it is evident that a magnetic inhomogeneity is causing 
triplet paring with parallel spins. Long range triplet currents 
should also appear in inhomogeneous magnetic junctions with 
(i-wave superconductors. 8 

Half-metals 9 are fully spin-polarized ferromagnets, i.e. 
they are metallic for one spin-direction and insulating for the 
other. Some half-metals include CrC»2, Lao.7Si"o.3Mn03, and 
Fe 3 04. 10 Keizer et a/.^were the first to report long-range su- 
percurrents through CrC<2. However, this group reported large 
variations in the magnitude of the critical current in the differ- 
ent samples. Moreover, some samples did not show any long- 
ranged supercurrent at all. Recent experiments have detected 
stable long range currents believed to be caused by a forma- 
tion of spin-triplet pairing. 12 All experiments report some sort 
of spin-flip scattering caused by a magnetic inhomogeneity is 
required to obtain the long-range currents. 

As there is only one spin direction at the Fermi level in 
half-metals, the current passing through will be completely 
spin-polarized. This has the potential for useful applications 
in low- temperature nanoelectronics. There have been several 
theoretical works on half-metallic Josephson junctionsfSKS 
but to our knowledge no-one has considered d-wave super- 
conductors in this context. From a more general perspective, 
the interplay between unconventional superconductivity, such 



as d-wave pairing, and half-metallicity has been studied in the 
context of magnetoresistance and spin-injection properties of 
cuprate/manganite hybrid structures 1922 . The study of a half- 
metallic d-wave Josephson junction therefore has relevance 
for these types of materials. 

In this paper, we will study how the supercurrent through 
unconventional half-metallic Josephson junctions depends on 
the properties of spin-active interfaces. We also study how the 
critical current depends on parameters like junction width and 
temperature for different levels of spin flip scattering. Com- 
paring the different superconducting symmetries, we find that 
a d^-wave superconductor gives a somewhat different result 
than s-wave and d a .2_ J ,2-superconductors. Junctions made of 
of half-metals and the latter two types of superconducting 
order show a non-monotonic temperature dependence of the 
critical current, while corresponding junctions involving d xy - 
superconductors do not show this behavior. Moreover, the 
junctions involving d^-superconductors may sustain a crit- 
ical current which is considerably larger than the two other 
superconducting symmetries. We also show that one can in- 
duce — 7T -transitions^ by switching the alignment of the 
magnetic interface fields from parallel to anti -parallel. It is 
also possible to tune the ground state superconducting phase 
by rotating the magnetic misalignment field at one interface 
relative to the other interface. These results are obtained nu- 
merically. We also provide a general analytical expression for 
the Andreev bound-state (ABS) energy spectrum which con- 
firms our numerical findings. 

The rest of the paper is organized as follows. In section 
II, we introduce our model and notation. Our results are pre- 
sented and discussed in section III, and we present our con- 
clusions in section IV. We use . ? . for 4 x 4 matrices, and ^^ 
for 2x2 matrices. Boldface notation is used for 3 -vectors. 



II. MODEL AND FORMALISM 

We use a modified BTK (Blonder- Tinkham-Klapwijk) the- 
ory^ w hich takes into account arbitrary pairing symmetries 
and a spin-active S/F interface. We consider an S/F/S junction 
of width L, see Fig.[T] The spin-active interfaces may induce 
long-range triplet currents. The propagation of quasiparticles 




FIG. 1 : Schematic illustration of our system. A half-metal of width 
L is sandwiched between two superconductors which are considered 
as reservoirs. The superconductors may have either an s-wave or d- 
wave symmetry. The interface regions are assumed to be spin-active 
due to e.g. magnetic disorder or misaligned local magnetic moments. 



is described by the Bogoliubov-de Gennes equations 

Hifj = £l[) 

where ip is the eigenstate with energy eigenvalue e, 
Hamiltonian of the system is given by 



H = 



Ho + V iA(0)a2 
-iA(0)*a 2 ~Ho-V* 



(1) 



The 



(2) 



where 

V = 



2m 



-E F \a -h z ®{x)®{L-x)o3 



V x - iV, 



V X + iVy V l 



v ) [S(x) + S(x - L)] 



(3) 



Here, O(x) and S(x) denote the Heaviside step-function and 
the delta-function, respectively. The cr-matrices are the Pauli 
matrices, m is the effective mass of the quasiparticles in 
both the superconductors and the ferromagnet and Ep is the 
Fermi energy. We assume equal Fermi energies in the differ- 
ent regions of the junction. The superconducting gap is de- 
noted by A(0) = A a (T)g(6)[e lipL e(-x) + e iVR 9(x - £)], 
where g(9) accounts for the superconducting pairing sym- 
metry, A (T) is the temperature dependent gap amplitude 
and <Pl(R) is the phase of the left(right) superconductor. We 
consider the usual BCS temperature dependence Aq(T) = 



Ao tanh 1 1.74y/T c /T — 1 j where T c is the superconducting 

critical temperature. The exchange energy is h z , and its di- 
rection is parallel with the z-axis. We will consider the limit 
where the ferromagnet becomes half-metallic, i.e. h z —} Ep. 
The barrier magnetic moment constitutes a spin-dependent 




FIG. 2: The barrier magnetic moment at the interface and its mis- 
alignment angles >]/ and <j>. The bulk magnetization in the ferromag- 
net is assumed to be aligned with the z-axis. 

potential V = (V x , V y , V z ), where the components are 

V x = — pVo cos * sin 4> 

V y = — pVo sin <3> sin cf) 

V z = - pV cos 4> (4) 

and V a = Vq + <rV z . a = ±1 for spin-up and spin-down. See 
Fig. E] for an illustration of the barrier magnetic moment and 
its misalignment angles "J and <\>. The non-magnetic barrier 
potential is Vo, while p is the ratio between the magnetic and 
non-magnetic potentials, i.e. 



P = \V\/V 



(5) 



We assume that the bulk magnetization of the ferromagnet is 
aligned with the z-axis in Fig. [2] For p ^ with <f> = we 
have spin-mixing and with <\> ^ we also have spin-flip. For 
details on spin-mixing and spin-flip, see e.g. RefPl 

Solving the BdG-equations yields the wavefunction in the 
different regions of our system. 25 J We can have have four 
different incoming quasiparticles, electron like quasiparticles 
(ELQ) with spin up and down, and hole-like quasiparticles 
(HLQ) with spin up and down. For an incident spin-up elec- 
tron in the left superconductor, the wave function is 

i> L (x)=[u L (6 + ),0,0,v L (6 + )e-^}e* k +* 



+ri[u L (9-) 1 Q,O,v L {0-)e-^}e 
r£[O,u L (0_),C«L(Me-^NO] 



r ' e 

\-rl 



e 



k+x 

-ik+x 



hr^K 



(e + y^,0,0,u L (9 + )}e' 






(6) 



For this particular process, the coefficients r|, r^,rl, r^ de- 
scribe normal reflection, normal reflection with spin-flip, 



novel Andreev reflection and usual Andreev reflection, re- 
spectively. We note that the momentum parallel to the inter- 
face is conserved for these processes. 

The corresponding wave function in the right superconduc- 
tor is 

M*) =tt[un(0+)e^, 0, 0, v R (9 + ) C -^]e^ x 

+4M#_KW^O,O, UR (0_)]e-^, (7) 

where the t's are the transmission coefficients, corresponding 
to the reflection processes described above. We have defined 
e^± = A(6±)/\A(6±)\ with 9+ = 9 and 6L = n - 8, and 
ip = ipn — ifL is the phase difference over the junction. The 
parameter £ accounts for singlet or triplet pairing in the su- 
perconductors. Here, we will consider only singlet pairing, 
i.e. £ = — 1. Previous works have considered tunneling in 
p-wave superconductor/ferromagnet structure J23E21 The co- 
herence factors are defined as usual 



u(9) = 
v{9) = 



\ 



s/^WW 



\ 



V^-IAWI 2 



(8) 



k± = \J2m(Ep ± \A 2 — |Al(6')| 2 ) cos 9 is the wavevector 
for ELQs (fc+) and HLQs (fc_) in the left superconductor and 
q± is the corresponding wavevectors in the right superconduc- 
tor. The wave function in the half-metal is 



i>HM{x) =(ee ik °* + /e- | *!( x - i >)[l, 0,0,0] 
+{e'e ik i x + fe- lk ^ x - L) )[Q, 1, 0, 0] 
+ (ge-< x + he< (x ~ L) )[0, 0, 1, 0] 
+ (g'e~ lk * x + h'e lk ^ x - L) )[Q, 0, 0, 1], 



(9) 



where k° h 



y / 2m(Ep + ah z ± e) cos 9. 
All scattering coefficients can be determined by matching 
wave functions at the interfaces. The boundary conditions are 

<Mo~) = VWo + ) 

d x [^ L {x) - ^hm{x)]\ x=0 = 2mVz,Vi(0), (10) 



for the left interface and 

d x [lpHM(x) ~ i>R{x)]\ x = L 



2mV R i> R (L), (11) 



for the right. We will later use the parameter Z = 2mVo/kp 
as a measure of interface transparency. Here, kp = \j2mEp 
is the Fermi wave vector. After the differentiation in the 
boundary conditions, we take the half-metallic limit kZ . — > 0. 



From the boundary conditions we obtain a system of linear 
equations which yields the scattering coefficients. 

With the scattering coefficients at han d, we can use a gen- 
eralized version of the Furusaki-Tsukudi^offiH formalism to 
calculate the Josephson current 



tt/2 
Un -tt/2 



(0) 



k+(uj n ) + k (uj n ) 



oi(w„) - a 2 (w„) 03 (w„) - 04 (w„) 



k+(u n ) 



k (w„) 



, (12) 



where uj n = (2n + l)ir//3 are fermionic Matsubara frequen- 
cies with n = 0, ±1, ±2, ..., and n n = y/u% + |A L (6»)| 2 . /3 
is the inverse temperature. k + (u> n ), k~(ui n ) and aj(w n ) are 
obtained from fc + , fc~ and a,; by analytically continuing e to 
iuj n . a,i with ? = 1, 2, 3, 4 are the ordinary Andreev reflection 
coefficients for incoming ELQ with spin-up or spin-down, and 
incoming HLQ with spin-up or spin-down, respectively. The 
summation over the Matsubara frequencies is performed nu- 
merically. 

If we instead neglect the contribution from the incoming 
quasiparticle, we get a homogeneous system of linear equa- 
tions 



Ax = 0, 



(13) 



where x = [r|, r\, r\, rf v t\, t\, tj v t\] and A is an 8 x 8 ma- 
trix obtained by expressing the scattering coefficients in the 
half-metal by the scattering coefficients in the left and right 
superconductor. By requiring a non-trivial solution of this 
system, det A = 0, we find the Andreev bound state energy 
spectrum, Ei. 

From the Andreev bound states we can find the Josephson 
current for a short junction £/£ -C 1, where £ is the super- 
donducting coherence length, in the ordinary wajEl 



^) = 2e E^/(^ 



(14) 



where e is the elementary charge and /(-Ei) denotes the 
Fermi-Dirac distribution function. We find the critical current 
from I c = max^ |/(y)|. 



HI. RESULTS 

For conventional s-wave pairing we use g(9) = 1. For d- 
wave we use g(9) = 008(26* — 2a) where a = and 7r/4 
correspond to d x 2_ y 2- and d^-pairing, respectively. We will 
use the superconducting gap Ao as a unit of energy. The Fermi 
energy is Ep = 1000A , the interface transparency is Z = 1 
and we use T = 0.2T C unless otherwise stated. We consider 
equal gap amplitude for the different pairing symmetries to 
ease the comparison. Unless otherwise stated, we assume that 
the magnetic misalignment angles <j> and \& on both interfaces 
are equal. 
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FIG. 3: Critical current as a function of junction width for p — 
(dashed lines) and p = 0.5 (solid lines) with misalignment angles 
ty = 7r/2 and <f> = tt/2. Blue, green, and red lines indicate $-, 
d x 2_ 2-, and d xy -wave pairing symmetry, respectively. 



Figure[3]shows the critical current for all three pairing sym- 
metries for two values of p with misalignment angles <p = tt/2 
and ^ = 7r/2. The first value, p = 0, corresponds to no 
spin-flip. As expected, there is no current in this case. How- 
ever, when spin-flip scattering is present, p = 0.5, we see 
that a long-range current appears. We find that junctions with 
s- and d x 2_ y 2 -symmetry behave very similarly, with d x 2_ y 2 
having a slightly lower current magnitude. The critical current 
is proportional to the gap amplitude multiplied by a weighting 
factor depending on the Andreev reflection coefficients (see 
Eq. ( fT2] i) averaged over the angles of incidence 9. As the s- 
wave symmetry is independent of 9, its average is larger than 



for the d, r 2_ 



The d T 



,,_. ii,^ ,, xy -symmetry carries a current approxi- 
mately three times larger than s- and d x 2_ y 2 for our parameter 
choice. 

Figurefflshows how the critical current depends on the mag- 
netic misalignment angle </> with "J = tt/2 for three junction 
widths kL = 10, 100, 1000 which corresponds to widths of 
approximately 2, 20 and 200 nm, respectively. Notice that we 
obtain the maximal critical current at a misalignment angle 
(j> > tt/2 when the misalignment at both interfaces is paral- 
lel to each other. If we rotate the misalignment 180°, which 
is equivalent to inverting the exchange field h z —> —h z , the 
maximum will appear at cf> < tt/2. As long as h z ^ 0, the sys- 
tem is not invariant under a spatial inversion. 34 Therefore, the 
direction of the exchange field relative to the misalignment 
field will determine whether the maximal current appears at 
misalignment angles smaller or larger than tt/2. 

If the magnetic moments at the interfaces are switched from 
parallel (p^ — pn) to anti-parallel alignment (p^ = — /Or), the 
current-phase relation changes sign as seen in Fig.B] This in- 
dicates a 0— 7r-transition. Hence, it is possible to induce 0— tt- 
transitions in half-metallic junctions by switching the relative 
direction between the two misalignment fields. Note that the 
other superconducting symmetries show the same behavior. 




FIG. 4: Critical current as a function of the misalignment angle (j>. 
Solid, dashed and dash-dotted lines correspond to junction widths 
kL — 10, 100, and 1000, respectively. Panels (a)-(c) show s-, 
d x 2_ y 2-, and d^-wave pairing symmetry, respectively. 




FIG. 5: Current-phase relation for a junction with d xy - 
superconductors on both sides. The black solid line shows paral- 
lel misalignment at both interfaces, while the red dashed line shows 
anti-parallel misalignment. 



We also note that for anti -parallel alignment, the maximal crit- 
ical current is obtained at a misalignment angle <j> — tt/2 for 
both interfaces. 

It is also possible to continuously change the ground state 
phase of the junction by rotating the misalignment field at one 
interface relative to the other. Fig. [6] shows the current-phase 
relation for a d xy -junction when the magnetic moment at the 
left interface points in the x-direction while at the right inter- 
face the magnetic moment is rotated in the xy-plane from be- 
ing parallel to the left interface, to pointing in the y-direction. 
Here, 4>l = 4>r = tt/2. We see that the current changes 
continuously with ^/j. Hence, we can obtain a finite current 
at zero superconducting phase difference by tuning the mag- 
netic misalignment fields. This is the same result as reported 



in RefP^for s-wave superconductors. We get the same results 
for s- and d x 2 _ y i -wave junctions. 

To gain further insight into the physical mechanisms under- 
lying our numerical results, we have solved Eq. ( fT3) analyt- 
ically. This yields a rather large and complicated expression 



for the ABS-energies, but it is nevertheless possible to infer 
how the interplay between the misalignment angles and su- 
perconducting phase difference is manifested. The resulting 
expression is 



£ = ±|A(0)| 




A 2 t cos (ip + x)± \JB\ + B 2 t cos (p + x) + B3T 2 cos 2 (p + \) 

C 



(15) 



where A4, Bi and C are large expressions that depend on junc- 
tion parameters like junction width kL, the barrier magnetic 
moment V, the angle of incidence 8 and the wave vectors 
kl h . t = ap 2 Vy sin 2 <f>L where a = ±1 for parallel and 
anti-parallel misalignment, respectively, and x — ^R — ^ l 
is the difference between the azimuthal angle of the right and 
left interface magnetic moments. As the Jospehson current de- 
pends on the ABS-energy differentiated with respect to ip, see 
Eq. (14) , much of the qualitative behavior can be explained 
with this expression. First, if we have no spin-flip, i.e <$> = 
and/or p = 0, then r = 0. Every term containing ip is mul- 
tiplied by a r-factor, which shows that no spin-flip equals no 
current. This corresponds to the result we see in Fig. [3] where 
the current vanishes without spin-flip. Second, we see that 
a change of the magnetic misalignment on the two interfaces 
from parallel, a = 1, to anti-parallel a = —1 is the same 
as shifting <p — > p + n, i.e. a — 7r-transition like the one 
we see in Fig. [5] Lastly, the %-phase shows that it is possible 
to continuously change the ground state phase of the junc- 
tion by tuning the interface magnetic moments as shown in 
Fig. [6] since x effectively renormalizes the superconducting 
phase difference: p — > p + x- 

The most significant difference between the three gap sym- 
metries is the temperature dependence of the critical cur- 
rent. The Josephson current in the case of d XJ ,-symmetry de- 
creases monotonically with increasing temperature, while the 
Josephson currents for the two other symmetries show a non- 
monotonic behavior. Figs. [7] demonstrate the temperature de- 
pendence of the critical current for the different gap symme- 
tries. 

This unusual temperature dependence can be explained by 
considering the proximity-induced density of states in the 
half-metal. The critical current is proportional to the energy 
integral of the den sity of states multiplied by the slope of the 
Fermi distribution!^ 



0.15 



dE N(E)d E f(E) 



(16) 



where N(E) is a proximity-induced density of states in the 
half-metal, which vanishes when superconductivity vanishes. 
As shown in Refs |2£I22] the proximity-induced density of 
states for a half-metallic Josephson junction with s-wave su- 
perconductors has a peak at finite energies. Spin-flip pro- 
cesses at the interfaces create long-ranged spin-triplet pairs 
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FIG. 6: Current-phase relation for a junction with d xy - 
superconductors on both sides. The misalignment field in the xy- 
plane is rotated from parallel to perpendicular, i.e. ^r changes from 
to 7t/2 in steps of 7r/8 (left to right). * L = and p L = pn = 0.5. 
Here, <J)l — <f ) R = t/2. 



that penetrate the interface and create a peak in the density of 
states in the half-metal due to the formation of bound-states 
(in the ballistic limit). This peak is the origin of the non- 
monotonic temperature behavior of the critical current which 
we also obtain, see Fig. |7(a)| Subgap bound-states may be 
thermally activated and contribute to the transport at a tem- 
perature T* of the order of the energies of the bound-states. 
Thus, when the temperature T is increased towards T* the 
critical current can increase, although one should also take 
into account the fact that the superconducting order parame- 
ter is suppressed with increasing temperature. For the case of 
s-wave and (^2 _ y 2 -wave superconducting gaps, the required 
temperature T* is of the order of the gap-amplitude. For a 
dzy-symmetric gap, the bound-states tend to be at lower ener- 
gies close to the Fermi surface 25 , yielding a lower T*, and 
hence a critical current which essentially is monotonically 
decreasing with temperature. We emphasize that all of this 
is contingent on the presence of spin-active tunneling barri- 
ers. Fig. 7(b) shows that d x i _ y 2 -pairing yields the same non- 



0.05 r 



0.04 



0.03 



0.02 



0.01 




0.4 0.6 

T/T c 

(a)s-wave, kL = 100 



monotonic behavior as s-wave pairing, while Fig. 7(c) shows 
that da; j, -pairing decays in the usual monotonic way. 

We expect the same behavior for other superconducting 
symmetries as well: If the bound-states are at energies close 
to the Fermi level, i.e. the density of states peaks at or very 
close to Ep, we expect the usual monotonic temperature de- 
pendence. However, if the bound-states are at higher energies, 
i.e. the density of states peaks at higher subgap energies, we 
expect a non-mono tonic temperature dependence. 



IV. CONCLUSIONS 
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(b)d x 2_ y 2-v/ave, kL = 100 
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(c)d xy -wave, kL = 100 

FIG. 7: Critical current as a function of temperature for all three su- 
perconducting symmetries at junction width kL — 100. Notice that 
the d^-symmetry decreases monotonically, while s- and d x 2_ y 2- 
symmetries have a non-monotonic behavior. The legends show the 
misalignment angle <fi in units of ir. 



In conclusion, we have analyzed the dc Josephson effect in 
an S/HM/S junction with spin-active interfaces. The possi- 
bility of spin-flip scattering at the interface creates equal spin 
triplet Cooper pairs in the half-metal, such that a long-range 
supercurrent may be sustained. We have considered the role of 
unconventional superconducting symmetries, in particular d- 
wave, and found that d-^-symmetry yields the largest critical 
current magnitude, while s- and d x i_ y i -symmetry show sim- 
ilar behavior. In addition, the temperature dependence of the 
critical current is qualitatively different in the d xy -wa\e case. 
Namely, while the current shows a non-monotonic behavior 
for s- and d x 2_ y 2 -symmetry, it decays monotonically in the 
dxy-wave case. This may be explained by considering the 
proximity-induced density of states in the half-metal. When 
the density of states in the half- metal has more weight near the 
Fermi level, as in the d x „-wave case due to the hybridization 
of the bound-states at the S/HM interfaces, the temperature 
dependence of the critical current is monotonically decaying 
as the temperature increases. On the other hand, if the density 
of states has a peak far from the Fermi level, we expect a non- 
monotonic temperature dependence. We have also discussed 
how it is possible to switch between a 0- and 7r -junction, both 
continuously and abruptly, by controlling the magnetic prop- 
erties of the interfaces. Finally, we have found a general an- 
alytical expression for the Andreev bound-state energy spec- 
trum which confirms our numerical results. 
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